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ABSTRACT 

We have conducted a long-term y-band photometric monitoring of M33 
on 95 nights during four observing seasons (2000 - 2003). A total number 
of 6418 lightcurves of bright objects in the range of 14 - 21 mag have been 
obtained. All measurements are publicly available. A total of 127 new vari- 
ables were detected, of which 28 are periodic. Ten previously known non- 
periodic variables were identified as periodic, 3 of which are Cepheids, and 
another previously known periodic variable was identified as an eclipsing bi- 
nary. Our derived periods range from 2.11 to almost 300 days. For 50 variables 
we have combined our observations with those of the DIRECT project, ob- 
taining lightcurves of up to 500 measurements, with a time-span of ^ 7 years. 
We have detected a few interesting variables, including a 99.3 day periodic 
variable with a 0.04 mag amplitude, at the position of SNR 19. 
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1 INTRODUCTION 



System atic searches for variable objects in M33 started with the seminal work of 



Hubble 



( 19261 ). who found 42 variables in this galaxy. He identified 35 of the variables as Cepheids, 



with which he establis 
conducted since then 


led the extragalactic nature of M33. A few additiona 


studies have been 


Hu 


bble&Sa,ndaJl953. 


van den Rerffh et al .111 975. 


Sa^ndaffe & Carlson 



m 



I 



Kinman et al 



19871 ) using photographic techniques. 
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The first CCD-based variability search in M3 3 was conducted as part of the DIRECT 



project (jKahiznv et al 



Stanek et al 



19981 ) with the goals of finding Cepheids and 



Eclipsing Binaries (hereafter EBs) that will yield a better distance estimate to o ur neighbor- 
ing ga laxies (M31 & M33). The first DIRECT observational campaign of M33 (jMacri et al 



was carried out with the Whipple Observatory 1.2-m telescope and with the 
Michigan-Dartmouth-MIT Observatory 1.3-m McGraw-Hill telescope. Three 10.8' x 10.8' 
fields were selected to cover the central part of the galaxy, labeled M33A, M33B and M33C, 
located North, South and South- West of the center, respectively. Observations during 42 
nights, from September 1996 to October 1997, revealed 544 variables, including 2 51 Cepheids, 



47 EB s and 62 unclassified periodic variables, in the M33A and M33B fields flMacri et al 
2001bl ). 

J [ II IQ 

The DIRECT second observational campaign (iMocheiska et al.ll2001al lbfl aimed to follow 
two detached EBs in M33A and M33B. Observations were conducted with the Kitt Peak 



National Observatory 2.1-m telescope, usin g 10.4' x 1[ 



coordinates of the M33A and M33B fields of 



Macri et al 



.4' FOVs centered on the central 



()2001br i. Two separate runs, each 



of 7 nights, were conducted in October and November 1999. A total of 892 new variables 
were detected in both fields, increasing the number of DIRECT variables to 1436. 

In this study we wished to continue the DIRECT thorough search for variable stars in 
M33, with an emphasis on long time span, so we could identify variables with longer time 
scales, and in particular stars with double periodicity. We therefore performed a photometric 
monitoring of the same three fields with the 1-m telescope at the Wise Observatory on 95 
nights from September 2000 until November 2003. 

Ideally, we would combine our data with the individual measurements of the DIRECT 
project and search the combined data for new variables. However, the DIRECT individual 
measurements were publicly available only for the variables identified in M33A & M33B. 
Therefore we could combine the data only for the stars already identified by DIRECT as 
variables. We did that in order to get a better coverage of the periodicity and to try and look 
for additional modulations. For all the other stars we performed periodicity and variability 
searches independent of the DIRECT data. Although the precision of our measurements 
were inferior to those of the DIRECT project, the long time span enabled us to detect a 
few unknown variables, some of which are quite interesting. Among those are an optically- 
periodic X-ray source and a periodic variable at a SNR position. 

Observations and data processing are described in § |21 Section El describes the temporal 
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Table 1. Fields central coordinates (J2000.0) 



Field name 


RA 


Doc 


direct 1 


01:34:05.1 


+30:43:43 


direct2 


01:34:00.0 


+30:34:04 


directs 


01:33:16.0 


+30:35:15 



analysis techniques and results. It also describes the method of combining our data with the 
DIRECT data and compares our results with the DIRECT. We present a few interesting 
variables in § 0] and state concluding remarks in § 



2 OBSERVATIONS AND DATA PROCESSING 



2.1 Observations 



We observed M33 with the 1-m telescope at the Wise Observatory from September 2000 until 
November 2003. Observations were carried out with a standard Cousins- Johnson V filter us- 
ing a Tektronix 1024x1024 back-illuminate d CCD, with a pixe l scale of 0.696 ± 0.002 "/pixel 
and an 11.88' x 11.88' overall field of view (jKaspi et al.lll999l l. Exposure time was 900 sec- 
onds. 

M33 central region was covered by three fields — directl, direct2 and directs, similar to 
the DIRECT three fields. Table ^ lists the fields central coordinates, and Fig. Q shows the 
field boundaries on a Palomar Quick-K survey image of M33^. 

During the three years of the project M33 was observed on 95 nights, with about 30 
nights per typical observing season (In the season of 2002 - 2003 M33 was observed only on 
3 nights). We made an effort to observe each of the three fields twice per observing night, 
although some technical and/or weather conditions did not allow all six exposures to be 
acquired on all nights. A total of 617 exposures were obtained and are listed on Table |21 
They include 286 exposures of directl, 177 of direct2 and 154 of directs. 

On the night of Aug. 31, 2003, at the l ast stage of our monitorin g of M33, a nova 
eruption was discovered in the directl field (jCaneshalingam &: Lill2003^ . To follow-up the 



nova evolution, we o 



2003 (jShporer et al. 



Dtained 111 exposures for the directl field between Sep. 7 and Sep. 22, 



20031 ). Those images are included in our present analysis. 



^ The compressed files of the Space Telescope Science Institute Quick-Survey of the northern sky are based on scans of plates 
obtained by the Palomar Observatory using the Oschin Schmidt Telescope. 
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Figure 1. Boundaries of the three fields over-plotted on a Palomar Quick-V survey image of M33. North is up and East is to 
the loft. Coordinates are in the J2000.0 system. 



Table 2. Mid-exposure JD of all 617 V exposures obtained in the course of this survey. 

direct 1 

1801.290 1896.308 2134.414 2181.451 2305.187 2891.514 2899.316 2900.515 
1801.301 1898.196 2134.508 2181.461 2305.198 2891.525 2899.327 2900.526 



direct2 



1801.324 
1802.297 



1842.532 
1842.543 



1927.283 
1931.200 



2134.542 
2134.553 



2166.489 
2166.500 



2257.203 
2257.214 



2301.286 2904.384 
2301.297 2905.345 



dircct3 



1801.337 1829.289 1925.290 2131.554 2151.492 2259.242 2305.232 2908.536 
1801.348 1829.300 1927.295 2134.457 2166.521 2259.253 2305.243 2914.344 



Mid-exposure JDs are given as JD - 2450000, for each observed field separately. Only a small 
sample of Table 1^ is presented here. The entire table is available in the MNRAS electronic 



issue. 
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2.2 Photometry 

Photometric processing was performed using IRAF packages and routines^. Biasing and flat 
field correction were done with the CCDPROC package with exposures taken nightly. 

For object identification and astrometry we created for each field a reference image 
by combining the best 20 im ages. A total of 6487 objects were identified using the IRAF 



daof ind task IStetson 



19871 ). After X,Y coordinates were transformed to the coordinate 



system of each of the images we applied aperture photometry with the IRAF phot task in 
a fixed position mode. Other photometric techniques we experimented with did not result 
in significantly different results. Aperture radius was set to 3.5 pixels (= 2.4 arcsec) for all 
images, as this value is the typical PSF FWHM. 

To remove systematic effects in our data we applied the newly developed SysRem algo- 



rithm ( Tamuz et al, 



20051 ) . which can remove systematic effects in large sets of photometric 



lightcurves. SysRem succeeded to decrease the stellar scatter of the brightest objects by up 

to about 50%. 

Astrometry was performed using the USNO A-2.0 star catalog (jMonet et al.lll998^ and 
the reference image. CCD X,Y coordinates were transformed to equatorial FK5 with second- 
order polynomials. Residuals RMS were < 0.45" and no single coordinate residual exceeded 
1.0". 

Inst rumental magnitude was transformed to real magnitude using the catalog of|] 



Macri et al 



(2001a|). This catalog contains V magnitudes of all objects observed by the DIRECT first 



M33 observational campaign. A linear transformation from instrumental to real magnitude 
was derived using objects which were successfully astrometrically matched with catalog ob- 
jects, allowing a maximum distance of 2". Fig. |21 shows the magnitude difference between 
Wise and DIRECT magnitude vs. Wise derived magnitude for those objects. The increased 
scatter in the magnitudes difference for faint stars results from the increased uncertainty in 
the magnitudes of those stars. 

The RMS of the difference between Wise and DIRECT magnitudes is 0.3 - 0.4 mag, vary- 



i ng slig htly from field to field. This scatter is simil ar to the scatter p resented by 



( 2001aL Figure 7, 2001b, Figure 6) and found by 



T^fifi fit al 



Mocheiska et al 



^ IRAF (Image Reduction and Analysis Facility) is distributed by the National Optical Astronomy Observatories (NOAO), 
which are operated by the Association of Universities for Research in Astronomy (AURA), Inc., under cooperative agreement 
with the National Science Foundation. 
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Figure 2. Difference between Wise and DIRECT magnitudes vs. Wise magnitude for astrometrically matched objects. 



Table 3. Object catalog 



Desig. 


R.A. 


Dec 


V 


a 




hh mm ss.ss 


dd mm ss.s 


[mag] 


[mag] 


W30144 


01 32 55.66 


30 33 02.3 


20.63 


0.15 


W31679 


01 32 55.68 


30 34 27.7 


18.79 


0.03 


W31624 


01 32 55.69 


30 35 34.8 


17.61 


0.04 


W30656 


01 32 55.73 


30 31 41.4 


19.65 


0.08 


W31717 


01 32 55.73 


30 34 30.1 


18.48 


0.03 


W30559 


01 32 55.73 


30 39 31.3 


20.08 


0.12 


W31082 


01 32 55.74 


30 39 22.4 


20.36 


0.17 


W31078 


01 32 55.77 


30 39 12.7 


19.71 


0.08 


W31484 


01 32 55.81 


30 37 13.0 


16.42 


0.01 


W30901 


01 32 55.83 


30 35 24.4 


19.80 


0.09 



Columns contain (1) designation, (2) R.A. (J2000.0), (3) 
Dec (J2000.0), (4) V magnitude and (5) magnitude RMS. 
Only a small sample of Table |21 is presented here. The 
entire table is available in the MNRAS electronic issue. 



2.3 Photometric Results and ON-Line Data 

Fig. 121 presents tlie RMS vs. tlie averaged V magnitude for all objects. Poisson and readout 
noises are dominant for objects fainter than mag ~16. For brighter objects the scatter levels 
off and systematic noise at the level of ~7 mmag becomes dominant. For stars brighter than 
~ 14 noise level rises again due to CCD saturation. 

Table El lists the position, V magnitude and RMS of all 6418 objects. Lightcurves, con- 
sisting of a total of 899,509 individual measurements are publicly available at the Wise 
Observatory FTP server at wise-ftp .tau. ac . il : /pub/shporer/m33/, which can be ac- 
cessed through anonymous FTP. Each lightcurve is given as a separate text file, named 'Ic' 
with a 5-digit internal number extension. The listing of Ic. 12121 is given in Table 01 




-0.6 
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Figure 3. RMS against real V magnitude for all 6418 objects. The decrease of the Poisson noise and read noise as objects 
become brighter, up to a magnitude of ~ 16, is the dominant feature of the diagram. For brighter objects a residual systematic 
noise of ~7 mmag is observed. The increased scatter for objects brighter than ~14th magnitude is due to CCD saturation. 



Table 4. Example of an on-line 

lightcurvc text file. 



1809.4479 


20.511 


0.146 


1809.4589 


20.591 


0.158 


1813.5006 


20.358 


0.131 


1819.2919 


20.490 


0.169 


1819.3029 


20.697 


0.184 


1819.3735 


20.522 


0.113 


1819.3844 


20.493 


0.120 


1820.2384 


20.437 


0.243 


1820.2494 


20.267 


0.170 


1820.3184 


20.595 


0.154 



Columns contain (1) mid- 
exposure JD - 2450000, (2) V 
magnitude and (3) magnitude 
uncertainty. Only part of the 
lightcurve is presented here, for 
guidance regarding form and 
content of all lightcurve text files. 

3 SEARCH FOR VARIABILITY 

Each Lightcurve was searched first for periodic modulation, and then for non-periodic vari- 
ability. For stars with available DIRECT data, we also applied periodic analysis to the 
combined data. 
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3.1 Periodicity Detection 

Periodicity search was applied to all lig htcurves with the Analysis of Variance (AoV) al- 
gorithm of ISchwarzenberg-Czernvl ()l989r i. For each trial frequency, u, data was folded with 
the corresponding period and then binned into 10 bins, using equal as possible number of 
points per bin. Two variances were calculated: Binned lightcurve variance, s^(z/), and the 
sum of bins internal variance, sl{i>). Periodogram value, 5'(z/), was taken as the ratio of those 
variances, sKu) / 32(1^) ■ 

For each lightcurve, we defined S as the value of the highest periodogram peak divided 
by periodogram average: 

_ maXyS{v) 
<S{.) >/ 

We consider the value of S as an indicator of the significance of the detection of a periodicity 
within the lightcurve. 

To estimate the significance of the periodicity detection we computed S for 100 random 
permutations of every lightcurve in each of the three fields, obtaining an S distribution 
consisting of ~ 2 x 10^ elements per field. We defined HsiS) to be the percentage of randomly 
permuted lightcurves with higher S values. Stars with Ti.s{S) smaller than 0.01% were 
flagged as periodic. This threshold gives an expectation of one false detection for every 10^ 
lightcurves, or, 0.64 for our entire sample. 

A total of 113 periodic variables were detected. The period uncertainty, AP, was defined 
as the FWHM of the periodogram peak. In order to derive the amplitude we applied an 
iterative running-median to the phased lightcurve and took half the magnitude difference 
between maximum and minimum brightness. 

We have classified 45 periodic variables as Cepheids by examining the period, ampli- 
tude and shape of all periodic variables. In particular, we examined the relative duration of 
increasing and decreasing brightness, in order to detect stars with an increasing phase sub- 
stantially shorter than the decreasing phase. Comparing with the publicly available DIRECT 
catalogs and the SIMBAD astronomical database we find that 8 out of the 45 Cepheids are 
new, and 37 are previously known Cepheids. The 8 new Cepheids are listed in Table El and 
plotted in Fig. |31 Five of the new Cepheids were not identified before as variables, and 
the other 3 were classified as variables and had no known period. The 5 new Cepheids are 
positioned in our direct3 field and in the DIRECT M33C fi eld. Variables of this field were 



not reported by the DIRECT project, although Table 9 of 



MacrietaL 



( 2001bl reporting 
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Table 5. List of 8 new Cepheids. Five of the new Cepheids were not identified before as 
variables and the other 3 were classified as variables and had no known period. 



Desig. 


R.A. Dec. V 


a Amp. Period 


AP Comments 




hh mm ss.ss dd mm ss.s [mag] 


[mag] [mag] [days] 


[days] 



New Variables 



W30837 


1 33 27.66 


30 34 24.0 


19.52 


0.12 


0.15 


14.837 


0.073 


1 


W31339 


1 32 57.55 


30 38 46.5 


19.36 


0.15 


0.16 


18.553 


0.096 


1 


W31573 


1 33 16.45 


30 36 58.2 


19.65 


0.15 


0.18 


18.553 


0.091 


1 


W31027 


1 33 30.15 


30 38 04.8 


20.45 


0.34 


0.35 


19.920 


0.100 


1 


W30348 


1 33 30.40 


30 35 55.8 


19.67 


0.36 


0.45 


22.08 


0.13 


1 



New Periodic Variables 



W22288 1 33 57.56 30 38 44.8 
W20128 1 34 15.47 30 31 06.7 
W21681 1 33 39.92 30 35 08.2 



18.04 0.07 0.08 26.74 0.16 
18.93 0.08 0.10 28.01 0.17 
18.31 0.12 0.15 51.28 0.91 



Columns contain (1) Designation, (2) J2000.0 R.A., (3) J2000.0 Dec, (4) mean V magnitude, 
(5) magnitude RMS, (6) amplitude, (7) period, (8) period uncertainty and (9) comments. 
Comments: 1; Position in DIRECT M33C field. 



variables in fields M33A and M33B) includes 15 M33C variables. (The 37 known Cepheids 
are included in Table ITTl which lists all variables identified here.) 



In order to detect EBs we applied the Eclipsing Binary Automated Solver (jTamuz. Mazeh &: North 



2006 



Mazeh. Tamuz fc North 



20061 hereafter EBAS) to all periodic variables that were not 



identified as Cepheids. EBAS has its own built in goodness-of-fit estimator, the alarm — A, 
which we used to identify the true EBs. The advantage of the alarm over the classical x^, as 
a goodness-of-fit estimator, results from accounting for measurements order, by considering 
maximal series of consecutive, mean-subtracted measurements with the same sign. We man- 
ually inspected all objects with A < 1.0, in order to remove spurious EB classifications. We 
identified 12 systems as EBs, of which 5 are new, presented in Fig. El and listed in Table IHl 
Those 5 include four new variables and one new classification of a known, unclassified pe- 
riodic object (W21568). Table [3 lists the EBAS solution parameters of all 12 EBs detected 
here. The 7 already known EBs are included in Table 

Periodic variables not identified as Cepheids or EBs were noted as unclassified periodic. 
Those include 56 objects, 26 of which were not known before as periodic variables, listed on 
Table IHl (The 30 known periodic variables are included in Table ITTl) The 26 new unclassified 
periodics include 19 new variables (not known before as variables) and 7 previously known 
as non-periodic variables. A sample of 8 new unclassified periodic variables is presented in 
Fig.ini 

In addition, we have determined an im proved period, of 175 .4 days, for W11984 where 
the previously known period was 458 days (jKinman et al .11 198 71 ). 
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W30837- 14.84 - 19.52 



W31339- 18.55 - 19.36 



19.3 
19-4 

(0 

g 19.5 



^ 19-61 * 



•: 
• • • 



19.7 



• • • 



0.2 0.4 0.6 0.8 

W31573- 18.55 - 19.65 



19.4 
(0 19.6 
^ 19.8 
20 



0.2 0.4 0.6 0.8 1 

W30348- 22.08 -19.67 



O) 19.5 
(0 

E 

> ^° 



20.5 



• •• 



0.2 0.4 0.6 0.8 

W20128- 28.01 -18.93 




0.4 0.6 

Phase 



19.2 



19.4 



19.6 




0.2 0.4 0.6 0.8 

W31027- 19.92 -20.45 



20 



20.5 



21 



« « 



• • V 



0.2 0.4 0.6 0.8 1 

W22288- 26.74 - 18.04 



17.9 



18.1 



18.2 




0.2 0.4 0.6 0.8 

W21681 -51.28- 18.31 




0.4 0.6 

Phase 



Figure 4. Lightcurves of the 8 new Cepheids detected here. Magnitude is plotted against phase. Title of each lightcurve is 
composed of (from left to right) designation, period, in days, and mean V magnitude. 



3.2 Non-Periodic Variability Detection 

All non-peri odic lightcurves wer e searched for variability using the alarm statistic, A, in- 



troduced by 



Tamuz et al. 



( 2006r i. Here, the alarm is used to detect variable lightcurves by 
estimating the goodness-of-fit of a constant function. The significance of the detection of the 
variability was estimated by a permutation test, as was done for the periodicity detection. 
For each single lightcurve, we generated 10^ random permutations and calculated their alarm 
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Table 6. List of 5 new eclipsing-binaries. Three of those 
were not identified before as variables, and the other one was 

idoiitifiod as an uiiclassifiod y)orioflic. 



Desig. 


R.A. 


Dec. 


V 




Comments 




hh mm ss.ss 


dd mm ss.s 


[mag] 


[mag] 




New Variables 


W30358 


1 32 57.31 


30 36 07.6 


19.04 


0.07 


1 


W31523 


1 33 29.88 


30 31 47.5 


17.52 


0.07 


1 


W30602 


1 33 35.31 


30 40 24.7 


20.15 


0.19 


1 


W20549 


1 33 58.67 


30 35 26.6 


16.73 


0.03 








New C'lassilical 


ion 






W21568 


1 33 54.80 


30 32 49.0 


18.52 


0.04 





Columns contain (1) Designation, (2) J2000.0 R.A., (3) 
J2000.0 Dec, (4) mean V magnitude, (5) magnitude RMS 
and (6) comments. 

Comments: 1: Position in DIRECT M33C field. 



Table 7. EBAS orbital solutions for the 12 eclipsing-binaries identified here. 



Desig. A V Period rt k Js x ecosw esinw Ap A, 

[mag] [days] 



W21612 


-0.09 


19.4681 


2.337324 


0.706 


0.2 


1.9 


0.79 


0.014 


-0.018 


0.995 


1.0000 






d=3.2c-02 


±8.5c-05 


±9.7c-02 


±1.2c+00 


±2.6e-|-00 


±1.3c-01 


±3.1e-02 


±4.4e-02 


±2.6c-02 


±2.3c-03 


W11486 


0.12 


18.359 


2.70811 


0.81 


0.93 


0.75 


0.800 


0.034 


-0.037 


0.99 


1.000 






±3.1c-02 


±l.le-04 


±1.0c-01 


±5.8c-01 


±2.7c-01 


±8.1c-02 


±2.4c-02 


±4.5e-02 


±2.1c-01 


±1.4c-02 


W11491 


-0.22 


18.980 


3.84684 


0.640 


1.44 


1.27 


0.64 


0.002 


0.08 


1.0000 


1.00 






±7.2e-02 


±3.2e-04 


±9.7e-02 


±8.0e-01 


±8.4e-01 


±l.le-01 


±4.7e-02 


±1.0e-01 


±3.9e-03 


±1.7e-01 


W10764 


0.55 


20.318 


4.432920 


0.731 


0.59 


1.43 


-0.09 


-0.060 


0.05 


0.016 


0.08 






±7.0e-02 


±7.1e-05 


±5.7e-02 


±5.6e-01 


±6.0e-01 


±1.3e-01 


±2.3e-02 


±1.0e-01 


±1.3e-02 


±1.6e-01 


W20942 


0.53 


20.739 


5.09566 


0.670 


1.06 


0.67 


-0.001 


-0.015 


-0.06 


1.0000 


0.34 






±4.9e-02 


±1.6c-04 


±6.8c-02 


±2.3c-01 


±1.9c-01 


±8.5e-02 


±2.5c-02 


±1.2c-01 


±9.2e-03 


±2.6e-01 


W30358 


0.16 


18.990 


5.48325 


0.634 


0.68 


1.36 


0.661 


0.047 


0.061 


0.000 


1.000 






±3.7e-02 


±3.3c-04 


±3.7c-02 


±4.4c-01 


±4.7c-01 


±5.2c-02 


±1.2c-02 


±5.6c-02 


±3.9c-02 


±1.6c-03 


W30051 


0.78 


18.5282 


6.624649 


0.6751 


2.64 


0.67 


-0.015 


-0.0645 


0.024 


0.07 


1.00 






±7.3e-03 


±l.le-05 


±9.7e-03 


±9.0e-01 


±1.2e-01 


±5.1e-02 


±5.4e-03 


±1.4e-02 


±2.0e-01 


±1.4e-01 


W30602 


-0.16 


19.980 


6.85114 


0.785 


0.93 


0.99 


0.41 


-0.017 


-0.106 


1.00 


0.00 






±6.6e-02 


±4.1e-04 


±6.1e-02 


±5.1e-01 


±1.6e-01 


±1.0e-01 


±1.7e-02 


±9.7e-02 


±2.8e-01 


±2.6e-01 


W31523 


-0.18 


17.473 


6.92557 


0.634 


1.71 


0.73 


0.784 


-0.0056 


-0.106 


0.005 


1.000 






±1.6c-02 


±1.2c-04 


±3.3e-02 


±4.0e-01 


±6.6c-02 


±1.7c-02 


±5.4c-03 


d=1.6c-02 


±3.4c-02 


±5.7e-02 


W22214 


-0.23 


18.582 


8.77401 


0.455 


3.8 


0.89 


0.32 


-0.031 


-0.001 


0.902 


0.99983 






±2.4c-02 


±4.1e-04 


±1.5c-02 


±1.3c+00 


±3.7e-01 


±1.7c-01 


±1.0c-02 


±3.9c-02 


±2.6c-02 


±1.5e-04 


W21568 


-0.20 


18.493 


9.7870 


0.549 


4.0 


0.050 


0.947 


0.002 


-0.067 


1.000 


0.13 






±2.9c-02 


±1.3c-03 


±3.9c-02 


±l.lc-|-00 


±7.2c-02 


±5.4c-02 


±3.1c-02 


±3.7c-02 


±3.3c-02 


±2.2c-01 


W20549 


0.27 


16.71 


24.4920 


0.652 


1.16 


1.343 


0.7938 


0.2464 


0.0600 


0.940 


0.36 






±l.le-l 


±3.9e-3 


±2.1e-2 


±l.le-l 


±5.2e-2 


±9.2e-3 


±6.7e-3 


±8.1e-3 


±6.1e-2 


±1.0e-l 



Columns contain (1) Designation, (2) alarm, (3) out-of-cclipsc V magnitude, (4) period, in days, (5) relative sum of radii, (6) 
radii ratio, (7) surface brightness ratio, (8) impact parameter, (9) eccentricity, e, multiplied by cosw and (10) by sinui, where w 
is the periastron longitude, (11) bolometric reflection coefficient for the primary binary component and (12) for the secondary 
component. 



value. For each real lightcurve we defined Ha{A) to be the percentage of randomly permuted 
lightcurves with higher A values. Lightcurvcs with 7Y^(A) smaller than 0.01% were flagged 
as non-periodic variables. 

A total of 177 non-periodic variables were detected. Those include 99 new variables, 
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Figure 5. Lightcurves of the 5 new EBs detected here. Magnitude is plotted against phase. EBAS Hghtcurve solution is 
presented by a solid line. Title of each lightcurve is composed of (from left to right) Designation, EBAS period, in days, and 
EBAS out-of-eclipse V magnitude. 



listed on Table El (The 78 previously known are included in Table ITTl ) A sample of eight 
new non-periodic variables is presented on Fig. [3 

Fig.lSlpresents a mean V magnitude histogram of all variable (periodic and non-periodic) 
objects detected here, indicating a completeness magnitude of about 19.5. For comparison, 
the corresponding histogram for all objects, where bin hight was reduced by a factor of 40, 
is also presented, in gray. 



3.3 Combined Lightcurves 



We c ombined the data pr esent ed in this paper with the available data for DIRECT variables 



from 



Macri et al. 



(|2nniah and 



Mocheiska et al 



Data from all three sources were 



combined to create < 7 years time span lightcurves, although with large gaps, of up to 



Wise Observatory M33 V -Band Variability Search 



Table 8. List of 26 new unclassified periodic variables. Twenty of those were not identified 
before as variables. The other 8 are known variables without a known period. 



Desig. 


R.A. 


Dec. 


V 


f7 


Amp. 


Period 


AP 


Comments 




hh mm ss.ss 


dd mm ss.s 


[mag] 


[mag] 


[mag] 


[days] 


[days] 




New Variables 


W10109 


1 34 15.57 


30 41 10.1 


18.61 


0.04 


0.04 


9.718 


0.021 




WoUo44 


i OO 2 i .6 I 


oO OO 51.5 


20.81 


0.24 


0.25 


11.211 


0.021 


1 


W30250 


1 33 25.58 


30 34 27.0 


20.32 


0.18 


0.19 


11.442 


0.047 


1 


W30917 


1 33 18.48 


30 35 34.3 


20.19 


0.18 


0.18 


12.107 


0.061 


1 


W30558 


1 33 31.65 


30 39 32.0 


20.38 


0.26 


0.27 


13.175 


0.048 


1 


W30569 


1 33 32.17 


30 39 46.2 


20.52 


0.29 


0.33 


14.948 


0.046 


1 


W20380 


1 34 16.05 


30 33 44.9 


16.79 


0.02 


0.03 


15.773 


0.063 




W12052 


1 34 00.90 


30 40 25.1 


18.25 


0.04 


0.05 


26.46 


0.34 




W12073 


1 33 45.15 


30 44 19.0 


18.02 


0.04 


0.04 


52.9 


1.2 




W31739 


1 33 12.34 


30 38 48.7 


17.54 


0.05 


0.05 


68.03 


0.80 


1 


W31230 


1 33 11.16 


30 34 21.9 


16.80 


0.03 


0.04 


99.3 


4.4 




W11579 


1 33 48.32 


30 42 41.4 


17.03 


0.03 


0.03 


104.2 


3.1 




W31005 


1 33 05.77 


30 37 20.3 


18.57 


0.09 


0.15 


105.3 


3.4 


1 


W31047 


1 33 15.76 


30 38 22.1 


17.96 


0.05 


0.07 


117.6 


3.8 


1 


W30803 


1 33 00.09 


30 33 41.4 


19.27 


0.10 


0.09 


140.8 


4.8 


1 


W21460 


1 33 50.96 


30 38 19.4 


16.97 


0.05 


0.08 


147.1 


5.0 




W30292 


1 32 59.34 


30 35 05.1 


17.91 


0.04 


0.04 


222 


13 


1 


W11898 


1 34 10.91 


30 41 41.6 


18.70 


0.05 


0.06 


270 


34 




W31416 


1 33 12.60 


30 32 52.5 


19.51 


0.20 


0.31 


278 


29 


1 


New Periodic Variables 


W30043 


1 33 36.60 


30 31 55.0 


19.36 


0.09 


0.09 


13.624 


0.068 




W21058 


1 33 41.34 


30 32 12.9 


18.24 


0.06 


0.06 


20.12 


0.16 




W21193 


1 33 58.33 


30 34 29.7 


18.04 


0.04 


0.04 


36.50 


0.29 




W21387 


1 34 05.42 


30 37 19.7 


19.16 


0.10 


0.13 


135.1 


6.5 




W20509 


1 33 55.58 


30 35 01.0 


17.39 


0.04 


0.06 


144.9 


7.0 




W11870 


1 33 52.39 


30 39 08.4 


16.55 


0.10 


0.15 


263 


10 




W10678 


1 33 58.93 


30 41 38.5 


17.11 


0.02 


0.04 


286 


26 





Columns contain (1) Designation, (2) J2000.0 R.A., (3) J2000.0 Dec, (4) mean V magnitude, 
(5) magnitude RMS, (6) period, (7) period uncertainty ,(8) amplitude and (9) comments. 
Comments: 1: Position in DIRECT M33C field. 

Table 9. List of 99 new non-periodic variables. 



Desig. 


R.A. 


Dec. 


V 


cr 


Comments 




hh mm ss.ss 


dd mm ss.s 


[mag] 


[mag] 




W31624 


1 32 55.69 


30 35 34.8 


17.61 


0.04 


1 


W31074 


1 32 56.73 


30 39 04.0 


19.22 


0.12 


1 


W31279 


1 32 57.86 


30 35 55.1 


19.13 


0.10 


1 


W30354 


1 32 58.19 


30 36 06.4 


18.80 


0.17 


1 


W31223 


1 32 59.58 


30 34 05.8 


19.23 


0.08 


1 


W30648 


1 32 59.70 


30 31 37.0 


19.84 


0.12 


1 


W31070 


1 32 59.74 


30 38 55.1 


18.06 


0.06 


1 


W30214 


1 33 00.94 


30 34 04.0 


19.19 


0.08 


1 


W30877 


1 33 01.02 


30 35 00.7 


18.79 


0.08 


1 


W30185 


1 33 01.70 


30 33 29.6 


19.55 


0.10 


1 



Columns contain (1) Designation, (2) J2000.0 R.A., (3) J2000.0 Dec, 
(4) mean V magnitude, (5) magnitude RMS and (6) Comments. 
Comments: 1: Position in DIRECT M33C field. 2: Position outside all 
three DIRECT fields. 

Only a small sample of Table |5]is presented here. The entire table is 
available in the MNRAS electronic issue. 
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Figure 6. Lightcurves of eight new unclassified periodics, out 
plotted against phase. Title of each lightcurve is composed of 
magnitude. 
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W30569- 14.95 -20.52 
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W11579- 104.20- 17.03 
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0.4 0.6 

Phase 

of the 26 new periodic variables identified here. Magnitude is 
(from left to right) designation, period, in days, and mean V 



500 V measurements each. Combining data from several telescopes requires fine calibration, 
removing possible zero point differences, when searching for variables in particular. 

To derive a periodogram we fitted the data from the three telescopes with a few harmonics 
for each frequency together with different zero points for the three data sources. Those zero 
points were fitted separately for each variable. Since there is a large scatter in the magnitude 



difference of the same stars in different telescopes (See Fig. Oland lMocheiska et al. 



( 2001a . 
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Figure 7. Lightcurves of eight new non periodic variables, out of the 99 new non-periodics detected here. Magnitude is plotted 
vs. JD - 2450000. Title of each lightcurve is composed of (from left to right) designation and mean Y magnitude. Data of only 
the first two observational seasons is presented on the above plots, to maJje the variability more apparent. 



Figure 7, 2001b, Figure 6)) taking the same zero point shift for all hghtcurves would not be 
accurate enough. 

Periodogram value for each frequency was taken to be the fitted amplitude divided by the 
goodness-of-fit parameter. We then divided the whole periodogram by a fitted polynomial, 
using a cubic smoothing spline fit. Number of harmonics used for each frequency in the fit 
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17 18 19 

V [mag] 



22 



Figure 8. Magnitude histogram of all variables detected here, showing a completeness magnitude of about 19.5. Corresponding 
(scaled) histogram for all objects is presented in gray. 



Akaike 



19741 ) ■ Fig. H shows 



was determined by using the Akaike Information Criterion (AIC, 
an example of this analysis for W10821. 

Lightcurves of 59 periodic objects were combined with data obtained here. For 9 of those 
stars results were unsatisfactory, mainly due to high noise level of faint objects. Results of 
the combining analysis for 50 objects are presented on Table ^1 together with the periods 
obtained previously by the two DIRECT project programs. 

Except for two objects, W10034 & W20923, period derived by the combining procedure 
is consistent with those derived by analyzing each data set separately. For W10034, analysis 
of the combined lightcurve was able to reveal the true period (107.53 days), which is twice 
the period derived previously. For W20923, we have detected a period of 69.44 days, while 
the DIRECT perio d is 56. 8 4. It is also interesting to note that a period of 69.50 days was 



already derived by 



Hubblel (jl926L Table II, variable No. 10) for this object. 



Table ^2 hsts all variables detected here along with designation of known variables. 



3.4 Comparison with the DIRECT results 

Many objects classified here as variables were not classified as such by the DIRECT project. 
Out of the total 290 variables detected here only 158 variable objects were astrometrically 
matched to DIRECT variables. The other 132 variables consist of 102 non-periodic and 30 
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Table 10. List of 50 combined lightcurves. 




DIRECT 




Combined 








Desig. 


■ 


PpirioH9 
J^cl 1UQ.Z 


Period 


AP 


Amp. 


ZPl 


ZP2 




[days] 


[days] 


[days] 


[days] 


[mag] 


[mag] 


[mag] 


W9 1 1 9 


2.3372 


2.3362 


9 '^'^797 
Z.OO 1 Z 1 


O. I a-yj^ 


n 1 s 


n of; 
u.uu 


n AO 


W1 1 A9.f, 


2.7081 


2.7104 


9 7ns 1 Q 
z. J uoiy 


u . yti-U'y: 


nil 


n 97 

u.z i 


n f\A 


W1 C\7(\A 
VV lU 1 yj^ 


4.4330 


4.4363 


4 43312 


2 Oe-04 


14 


U.U ( 


n SI 

U.Ol 


VV zuy^z 


5.0950 


5.0920 




1 6e-04 


n 97 

u.z / 


17 


71 

U. / 1 


W91 (\f\7 

VV Z J_UU / 


_ 


5.298 


o.zuio 


1 4e-03 


n f\A 

U.U4: 


-u.uo 


n 97 

U.z i 


VV ±UOU'4: 


8.14 


8.14 


O. -LOU / 


A (T^ 
4:.Oti-UO 


22 


n 1 Q 
u. ly 


u.uo 


W91 ^fiS 

VV Z J_OlJO 


_ 


4.869 


y . 1 oy oz 


Q 7r> OA 


n c\A 

U.U4: 


n DA 

-U.U4: 


n 1 s 

U. -Lo 


W1 C\OAA 
VV lUZ^^ 


11.95 


11.97 


1 1 Qf^SQ 
1 ±.yijoy 




ft "iA 
U.04: 


u. -Ly 


n 9Q 

-u.zy 


W1 1 9(^fi 

VV 1 J-ZUo 


12.35 


12.36 


IZ.OtJO 


5 le-02 


n 1 7 


nil 

U. -L 1 


n OS 

u.uo 


W1 1 f^79 
VV 1 -l-O i z 


12.92 


12.93 


19 Q1 QQ 
iz.y j_yy 


y . oti-uo 


n 9fi 

u.zu 


-u.uu 


u.uo 


wi ni 9fi 

VV ±UlZO 
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13.04 


1 (\(\'^A 


4:.ot;-UO 


n 9S 

u. zo 


u.uo 


n '?s 

u.oo 


VV zuyou 


13.31 


13.31 


i-O.O L:J\J 


O.Oc-UO 


U.Ou 


nil 

-U. 1-L 


u.uu 


W99'^4'^ 


13.56 


13.56 






n 9fi 

u.zu 


U. -LO 


-u.uy 


VV zuo t o 


14.35 


14.35 


1 A '\A7 


1 6e-02 


u.oo 


n 90 

u.zu 


n f^l 

-U.Ol 


WI 1 9Q1 
VV 1 J-Zy ± 


14.59 


14.61 




7 9r> 0*^ 


21 


nil 

-U. 1 -L 


-u.uy 


WI nSQ9 
VV luoyz 


14.92 


14.93 


1 A Q9^ 


1 3e-02 


9fi 
u.zu 


n on 
u.uu 


n 9Q 

-u.zy 


W91 A9S 
VV Z J_^ZO 


15.08 


15.08 




Q Ho 09 


n 98 

u.zo 


-0 12 


n A7 


VV zuyoo 


15.78 


15.78 




1 Oe-02 


u.oo 


n 07 

-U.U / 


n fii 

-U.Ul 


W991 

VV ZZ lOU 


16.28 


16.28 


1 (\ 97*^ 


6 le-02 


n "iA 


-0 14 


-0 21 


VV zuouo 


17.48 


17.48 




4 7e-01 


n 9Q 

u. zy 


n 1 Q 
u. ly 


-u.uo 


W11733 


17.83 


17.83 


17.841 


3.3e-02 


0.39 


-0.14 


0.27 


WI n9fi7 

VV ±uzu / 


18.61 


18.60 




7 Qr> fT^ 
/ . ye-uo 


u.ou 


n DA 


u.uu 


W1 1 704 
VV ± -1- 1 y± 


18.81 


18.81 


18.797 


1 7e-02 


0.50 


-0.04 


0.36 


WI nS91 

VV lUoZl 


20.19 


20.18 


20 161 


3 2e-02 


n 9fi 
u.zu 


n AS 


n 47 


W91 O'^S 

VV Z J-UiJO 


- 


- 


20.161 


1. le-02 


0.07 


-0.03 


0.11 


VV OUolJ- 


20.18 


20.18 


20 182 


1 8e-02 


24 


-U. lU 


n (19 

u.uz 


wi nn4fi 

VV iUUtU 


21.91 


21.93 


21.954 


1.9e-02 


0.40 


0.02 


0.22 


VV OJ-UlU 


22.67 


22.67 


99 fi^n 
zz . uou 


2 le-02 


u.ou 


n OA 


n (11 

U.Ul 


VV ZUUfi 


23.31 


23.31 


23.310 


1.6e-02 


0.45 


0.03 


0.27 


W9ni 'iS 

VV ZUlOO 


24.88 


24.88 




fi 7f> fT^ 

u. / e-uo 


u.oo 


(1 1 (1 

-U. lU 


n 1 Q 
u. ly 


WI DSO'^ 
VV ±uouo 


26.39 


26.37 


26.316 


2.7e-02 


25 


0.20 


-0.23 


WI riA'il 

VV lU^O-L 


26.48 


26.48 


9fi AQn 


3 2e-02 


n 9Q 
u.zy 


(1 (IQ 

u.uy 


n 'if) 
u.ou 


W91 97^1 

VV Z J_Z i u 


26.50 


26.50 


ZU. UZO 


4 4e-02 


u.ou 


(lis 

U. -LO 


U. lO 


WI ni 4*^ 
VV ±ui^o 


26.57 


26.58 


zu. oyu 


1 9e-02 


u.ou 


u.uu 


U. lO 


VV ±J-OOZ 


27.28 


27.30 


27 248 


1 8e-02 


n A'i 


-u.uo 


n (ifi 
u.uu 


VV lUUOU 


27.97 


28.00 


28 Oil 


3 2e-02 


24 


U. -LO 


-U. lO 


W'^nQ9fi 
VV ouyzu 


30.29 


30.29 


Qfi 91^7 


3 4e-02 


u.oo 


n c\A 


n 9S 

-u.zo 


W10773 


30.50 


30.52 


30.488 


5. le-02 


0.24 


0.16 


-0.05 


W10367 


31.25 


31.26 


31.348 


2. Oe-02 


0.46 


-0.06 


-0.03 


W10629 


33.95 


33.98 


34.014 


2.6e-02 


0.23 


0.13 


-0.05 


W11386 


35.86 


35.87 


35.84 


1.5e-01 


0.40 


-0.32 


0.17 


W21290 


36.09 


36.09 


36.036 


6. le-02 


0.51 


-0.28 


0.49 


W21743 


36.35 


36.35 


36.496 


4.8e-02 


0.19 


0.02 


-0.12 


W20627 


37.34 


37.34 


37.453 


2.9e-02 


0.47 


0.06 


-0.01 


W20702 


55.91 


55.91 


55.10 


1.5e-01 


0.33 


-0.01 


-0.17 


W10274 


55.99 


55.79 


55.25 


2.3e-01 


0.12 


0.00 


0.18 


W10466 


57.45 


57.45 


57.47 


1.2e-01 


0.16 


-0.06 


-0.18 


W20923 


56.84 


56.84 


69.44 


1.9e-01 


0.18 


-0.06 


0.20 


W20280 


73.94 


73.94 


74.63 


7.5e-01 


0.29 


-0.27 


0.46 


W10034 


55.0 


54.930 


107.53 


6.3e-01 


0.10 


-0.06 


0.31 



Columns contain (1) W ise Designation, (2) period of 'Ma ori et alJ I2OOIIJ1 . (3) period of 
iMocheiska et all l2001all3) . (4) period o f the combined lightcurves, (5) Period uncertainty, (6) am - 
plitude of the fit, (7) zero point shift of lMocheiska et al. 12001a b) relative to lMacri et al.l i2001hl) . 
and (8) zero point shift of the Wise data relative tO iMacri et al.i t2001bl) . 
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Figure 9. Combined lightcurve analysis of W10821 and its DIRECT counterpart — D33J013343. 4+304356.5. T wo harmonics 
wer e fitted to 462 measurements (of a time span of 2 615 days). Top panel sho ws data of the DIRECT 1st campaign iMacri et alJ 
l20blM . in circles, of the DIRECT 2nd campaign <Mocheiska et aljl20oTah . in crosses, and data obtained here, in triangles. 
Middle panel presents the periodogram where the 2nd strongest frequency is consistent with one half the strongest frequency 
value, respectively. Bottom panel shows the combined phase- folded lightcurve, where the fitted zero point of each data set 
was subtracted. Fitted zero points are 19.43, 19.91 and 19.90 mag for the DIRECT 1st campaign iMacri et alJl2001bl marked 
by circles), DIRECT 2nd campaign iMochciska et al. 2001a, marked by crosses) and the Wise data (marked by triangles), 
respectively. The black solid line is just an harmonic fit, derived by the lightcurves combining procedure. 



periodic variables. Of those 132, 5 were alread y id entified as varia 



flHubble &: Sandage 



1953, 



Kinman et al 



1987 and 



Shemmer et al. 



jles b y previous studies 
200ol ). giving a total of 



127 new variables, consisting of 99 non-periodic and 28 periodic. In addition, we obtained 
an improved variability type for 11 of the 158 astrometrically matched variables. 

The M33 area mon itored here includes thr ee fields, while the DIRECT catalogs of 

list variables from only two fields, M33A 



Macri et al, 



l20Qiy) and 



Mochejska et al 




and M33B. Those fields are also slightly smaller than the corresponding directl and di- 
rect2 fields observed here. Considering only the 221 variables we find here inside the FOV 
of DIRECT M33A and M33B fields, 148 of them are present in the variables catalogs of 
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Figure 10. Combined liglitcurve of W10034 witli its DIRECT counterpart — D33J013420.8+303943.0, revealing the true 
period of 107.53 days, twice tlien period derived previously. Only by analyzing the combined lightcurve, consisting of 490 
measurements taken in a time span of 2608 days, the true period was derived. The different markings for each data set are the 
same as in Fig. |3 bottom panel. 



Table 11. List of all 290 variables detected in this survey. Within each variability type objects are sorted by 
increasing period and the non-periodic variables are sorted by increasing R.A. For previously known variables their 
known designation is given. 



Desie 



R.A. 

hh mm ss.ss 



Wise 

Dec. 
dd mm ss.s 



Type 



Period 

[days] 



Designation 



Known 

Type Period Comments 
[days] 



W21612 
W11486 



01 33 56.26 
01 33 54.32 



+30 33 58.8 E 
+30 40 28.9 E 



2.337324 D33J013356.2+303358.6 E 2.34 

2.70811 D33J013354.3+304029.6 E 2.71 



W11572 01 33 45.93 +30 42 30.3 C 12.9199 D33J013346.0+304231.9 C 12.93 10 

W10128 01 34 02.76 +30 41 44.4 C 13.0634 D33J013402.8+304145.7 C 13.04 10 



W21162 01 33 60.96 +30 33 54.8 P 2.1142 D33J013359.8+303354.9 P 2.12 

W21067 01 33 41.63 +30 32 20.5 P 5.2618 D33J013341.6+303220.3 P 5.30 



10 



W31624 
W31074 



01 32 55.69 
01 32 56.73 



+30 35 34.8 
+30 39 04.0 



N 
N 



1,4 
1,4 



Key to variability types: E — Eclipsing binary, C — Cepheid, P — Unclassified Periodic, N — Non-periodic variable 
(DIRECT miscellaneous type). 

Comments: 1: New variable; 2: New classification of a known variable; 3: Position in DIRECT M33C, reported in 
Table 9 of Maori et al. (2001b); 4: Position in DIRECT M33C (new variable); 5: Position outside all DIRECT fields 
(new variable); 6: Object detected as periodic variable by Shemmer et al. (2000); 7: Astrometric matching between 
Wise and DIRECT objects can not be determined accurately due to the crowded region; 8: classified as a Long 
Period Variable (LPV) by Kinman, Mould & Wood (1987); 9: Object detected as variable by Hubble & Sandage 
(1953); 10: Period derived from the combined lightcurve. 

Only a small sample of Table [TTl is presented here, including two objects from each variability type. The entire table 
is available in the MNRAS electronic issue. 



Macri et al.l ((2001b|) or lMocheiska et al.l (|2001al jbl). The other 73 variable objects include 63 
non-periodic and 10 periodic variables. 

Table ^1 lists the number of variables detected here, of each variability type, and of 
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Table 12. Summary of the detections by the Wise Observatory 



Var. Type 


All Fields 


M33A S 


z M33B 


M33C 


Out 


Previous 


New Class. 


New Var. 






Wise 


DIR 


Wise 


DIR 










Cepheid 


45 


34 


34 


11 


6 






3" 


5 


EB 


12 


8 


7 


4 


1 






1*- 


4 


Unclass. Per. 


56 


40 


34 


16 


1 




2 


ya 


19 


Non Per. 


177 


136 


73 


38 


2 


3 


3 




99 


Total 


290 


218 


148 


69 


10 


3 


5 


11 


127 



Columns contain: (1) variability type and (2) number of detected variables. (3) Variables detected here 
inside the DIRECT M33A and M33B fields and (4) of those, number of variables detected by the DIRECT, 
including objects for which a new variability class was obtained here. (5) Variables dete cted here inside the 
DIRECT M33C field and (6) of those, number of variables detected by the DIRECT iMacri et al.ll2001bl 
Table 9). Only variables in M33C and not in overlapping regions are listed on columns 5 and 6. (7) Variables 
detected outsi de the DIRECT FOVs. (8) Variables detected here ^ detected previously by studies other than 
the DIRECT. <Hubble fc Sandagelll953llKinman et al.lll987l and lShemmer et alJbOOOt) . (9) Newly classified 
variables and (10) new variables. 

Notes: (a) classified by DIRECT as non-periodic, (b) classified by DIRECT as an unclassified periodic. 



those, number of variables detected by the DIRECT project. Table [T^ also refers to variables 
positioned in the DIRECT M33C field since although th is field's variables we re not reported 
by the DIRECT, 15 of them are included in Table 9 of iMacri et al.l ()2001bl ). Also given in 
Table IT^ are number of variables detected outside all three DIRECT fields, variables detected 
by other, previous studies, number of newly classified variables and of new variables. 



4 INTERESTING VARIABLES 

4.1 W21067 — An Optical Periodic X-Ray Source 

The combined data of W21067 (1:33:41.63, +30:32:20.5), one of the brightest objects in 
the sample, yielded a period of 5.26 days with an amplitude of 0.04 fse e Fig. ITTjl. Aii X-ray 
source at a distance of 1.4" from our position of W21067 was detected bv lPietsch et al.l ()2004 , 
Table 3, source 194) using the XMM-Newton observatory. They have measured an 0.2-4.5 
KeV flux of (1.3900 ± 0.0097) x 10~^^ erg/sec/cm^ and also suggested that the positional 
correlation between an optically periodic variable and an X-ray source makes this object an 
XRB candidate. 

To follow this suggestion Alceste Bonanos of the CfA obtained for us a multi-order 
spectrum of W21067 with the Echellette Spectrograph and Imager (ESI) at the Keck II 
10-m telescope in the Echelle mode. Fig. IT^ presents the joined spectrum of all 10 orders 
(orders 6 to 15), with a complete spectral coverage from 3900-10000 A. Velocity dispersion 
is about 11.4 km/sec/pixel in all orders. 

By correlating W21067's spectrum with the ELODIE library of template spectra we were 
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Phase 

Figure 11. Combined ligh tcurve of W21067 fo lded on the period detected by the combining proc edure, of 5.2618 days. D ata 
of DIRECT 1st campaign <Macri et aT]|2001hD is marked in circles, of DIRECT 2nd campaign <Mocheiska et alj|200l'bh in 

crosses and data obtained here in triangles. The solid line is a single harmonic fit calculated by the combining procedure (See 
section YA.'il . 




7000 

Wavelength [Angstrom] 

Figure 12. Joined spectrum from all 10 orders of W21067, taken by the Keck II ESI. Normalized flux is plotted against 
wavelength. Spectrum was smoothed using a running mean of 5 adjacent pixels. 



able to determine it is a K star, with a heliocentric radial velocity of +16.4 ± 0.6 km/s. The 
derive d velocity is significantly dif ferent from the systemic radial velocity of M33, of —179 ±3 
km/s ( de Vaucouleurs et al.l 1991 ). The most plausible conclusion therefore is that W21067 
is not an M33 member but a Galactic foreground object. 



22 Shporer & Mazeh 

The many absorption lines and absence of significant emission features in W 21067 spec- 



trurn rules out the possibility of a compact binary companion as the X-rays source (jBradt fc McClintock 



19831 ). Therefore we propose chromospheric activity as the source of the optical periodicity 



and X-ray radiation. 

The periodic modulations of chromospherically active stars are known to show long-term 



amplitude and shape variations, resultin 



activity cycles (|Guinan &: Gimenez 



g from their dark regions evoluti on during solar-like 



1992 



Olah et al. 



200G 



Rodon 3l992r) . This might be the 



cause for the somewhat large scatter in W21067 combined folded lightcurve, consisting of 
data from a period of ~7 years. 

Fill-in cores of the Ca II H & K absorption lines, at 3968.5 A and 3933.7 A, respectively. 



are a classical spectral feature used to identify chromospheric activity (iFekel fc Balachandran 



Strassmeier et a) 



19931 ). Fig. Cni shows W21067 Ca II H & K lines together with an 



active star and two inactive stars of a similar spectral type. Despite the decreased S/N 
close to the spectrum's blue end caused by the lower CCD QE and increased atmospheric 
extinction at these wavelengths, fill-in cores of W21067 Ca II H & K lines can be noticed. 



The color-temperature color-index relation (jAUen 



19731 ) for the DIRECT B -V = 
0.79 ± 0.04 color index yields a color temperature of Tc = 5250 ± 150 °K. For Galactic 
object in the direction of M33 extinction and reddening can be neglected, and therefore we 
can assume the ten aperature es timate is valid. The effective temperature bolometric correc- 
tion (BC) relation ( AUeJ 1973 ) yields BC = 0.18 ± 0.03 mag, which results in an X-ray to 
V flux ratio of logj fx I f v) = — 2.72 ± - 04, co nsistent with the mean value of —2.9 ± 0.8 



(median —2.8), given by 



Padmakar et al. 



( 20001 ) for 202 active binaries. 



4.2 W31230 — An Optical Periodic at SNR Position 



Gordon et al 



( 19981) and 



W31230 (1:33:11.17, 30:34:21.9) is located 0.41" from SNR 19 of 
was detected here as an unclassified periodic variable. Considering the number of SNRs 
and periodic variables detecte d here in the direct3 field, the probability of a random po- 



sitional correlation is 0.0009. 



Sholukhova et al 



Calzetti et al, 



( 1995[ ) list this object as an Ha emitter and 
( 1999( ) detected an Ha emitter at a distance of 2.2" from W31230, for 
which they measured an Ha equivalent width of 19 A. W31230 is located at the DIRECT 
M33C field and was not reported as variable by them. Our analysis detected a period of 
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Figure 13. From top to bottom: Spectrum of W21067 at the Ca II H & K spectral region. Spectrum of a chromospherically 
active star, HD 201091, a BY Dra K5 V star Spectra of two inac tive stars for comparison (HD 4628 K2 V, HD 16160 K3 V). 
Spectra of comparison stars were taken from lMontes et al.l |^93). 
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Figure 14. Top panel presents a Lomb-Scargle power spectrum of W31230 (normalized to unit mean) and the bottom panel 
shows the phased lightcurve, with an amplitude of 0.04 mag. 



99.3 ± 4.4 days and an amplitude of 0.04 mag. Fig. El presents a Lomb-Scargle (jScargle 
19821 1 power spectrum and a phased lightcurve of our data for this star. 
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4.3 Optical Variables at Wolf-Rayet Positions 

WR Variability can originate fro m random winds variat ions, rotation, spots evolution and 



pulsation, radial and non-radial (|Moffat fc Shara 



198^. A persistent periodic variability 



might indicate a WR in a binary system ^ yhere variability is ind uced by either geometric 



eclipses, wind eclipses or proximity effects (jMarchenko et al 



By correlating our variables against the SIMBAD database we found two objects which 
are located at positions of WR stars and a third object located at the position of a WR 
candidate. All three astrometric matches are within 1.1". 



W20549 (1:33:58.67, 30:35 :26.6), a spectroscop icallv confirmed M33 WR (jMassev fc Johnson 



1998L WR115) of type WNE (jMassev et al.lll99,5h . was identified here as a new EB, with 



a period of 24.4920 days and orbital eccentricity of 0.25 (see Fig. El and Tables IHl & IZ|). 
W31347 (1:33:27.25, 30 :39:09.7), a spectroscopically confirmed M33 WR of type WN9 



fMassev .Johnson 



1998L WR39), was detected as a new non-periodic variable ( see Fig. [3 



Massev et al. 



andT ableH. W10975 (1:34:20.30, 30:45:45.4) was identified as a WR candidate by 
()l987l Table 2, object 29) and was detected here as a new non-periodic variable (see Fig. [7| 
and Table El). 



4.4 W31284 - M33 Variable C 



W31284 (1:33:35.13 30:3 6:00.8) is a known LBV, named M33 Variable C, discovered by 



Hubble fc Sandagd ()l953L Fig. 6) by using data from 1921 to 1953. Sin ce then, this star was 



inclu ded in a few other lon g-term surveys: 
1972, 



R.osino fc Bianchinil fll973L Fig. 8 ), from 1960- 

( 1999[ Fig. 2) from 



Kurtev et al. 



Kinman et all ()1987L Fig. 5) from 1982-1985 and 
1982-1990. 

In both DIRECT campaigns and in this work M33 Var. C was identified as a non- 
periodic variable. Lightcurve containing data of the three telescopes, with a time span of 
7.16 years, is presented on Fig. ^1 Although there might be small zero-point shifts be- 
tween telescopes, a decrease of more than 1 magn itude, followed by a similar increase is 



evident. This lightcurve and the one presented by 



Kinman 



variations with a m u ch sh orter time-scale than observed by 



19871). show brigh tness 



Hubble Sandagel ()l 9531) and 



Rosino fc Bianchinil (|l973l ). Furthermore, Fig. 1151 shows that the star stayed at minimum 



brightness for a relatively short time before brightening again, a behavior which was not ob- 
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Figure 15. M33 Var. C lightcurve from DIRECT and this work, y-band data of the DIRECT first campaign is presented in 
circles, of the DIRECT second campaign in squares and from this work in triangles. DIRECT first campaign /-band data is 
marked by crosses. V-haad lightcurve shows a decrease of about 1 magnitude in the star's luminosity followed immediately 
by a similar increase. Variability time-scale is shorter than observed by Hubble & Sandage (1953) and Rosino and Bianchini 
(1973). 

served previously. This suggests that M33 Var. C brightness variations became more rapid 
since the 1980s. 




4.5 W11324 - LBV Periodic Microvariability Confirmation 



Masse V et aL 



W11324 (1:3 4:06.59, 30:41:46.8) is known a s B416 (jHumphrevs fc Sandagelll980f ). It was 
identified bv Shemmer. Leibowitz fc Szkodvi (|2000l hereafter SLS) as an LBV, confirming 



(19960 LBV candidate classification. SLS also report a periodic microvariability 



of 8.26 days with a V amplitude of 27 mmag, a V magnitude of 16.7±0.1 and an Ha EW 
of 106±2 A, c o nsiste nt with a V magnitude of 16.76 and Ha EW of 109.1 A, measured by 



Calzetti et al 



( 199,51 ). W11324 was not classified by DIRECT as a variable. Here, W11324 
was detected as an unclassified periodic, with a period of 8.251±0.017 days. The consistency 
of our period and SLS period, derived from observations taken more than 10 years apart, 
confirms the periodic nature of this object. 



4.6 W10913 - Hubble's Variable 19 in M33 



W10913 (1:33:57.02, 30:45:11. 5) was identified as variable bv lHubblel (|l926[ V19) who classi 



fied it clS Qb 54.7 days Cepheid. 



Macri et al 



( 20011 ) have shown that this object has undergone 



a dramatic decrease in its brightness variation, from a B amplitude of 0.55 mag in Hubble's 
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lightcurve to no detectable variation in the DIRECT data. In our data, this object is also 
non- variable and has a V RMS of 0.03 mag, similar to the scatter of other non-variable stars 
of the same brightness. 



5 CONCLUSIONS 



We presented here the results of a long-term monitoring of the bright stars in M33, complete 
to ~ 19.5 mag. We discovered 8 new Cepheids, 5 binaries, 26 unclassified periodic variables 
with period range up to almost 300 days, and 99 non-periodic variables. Combined with the 
publicly available data of the DIRECT project, our data covers more than 5 years of stellar 
variability. One of these variables is the famous M33 Var. C, which displays more than 1 mag 
modulation with a timescale of 2500 days. The variability found here is different from that 
previously reported and therefore further observations of this star will be most interesting. 
One intriguing result of this study is the findings abo ut W31230, one of th e brightest 



Cordon et a1 



( 19981 ) and shows 



stars in M33, which is found at the location of SNR 19 of 
Ha emission. It displays a small periodic modulation of 0.04 mag at a period of about 100 
days. This is si milar to one of the two periodic modulations found for the famous Calactic 



source SS 433 (iMargon 



Kirshner fc Chevalier 



1984^ . which is also located at the center of Calactic SNR (e.g.. 



19801 ) and has a prominent emission. The 165 d periodicity of SS 



433 is associated with the precessing disc around the compact object, which is probably 
the remnant of th e SN explosion. SS 433 displays anot her periodicitv of 13 davs (e.g.. 



Mazeh et al, 



19871 ). associated with the binary period (e.g.. lCherepashchuklll98ll ). We could 



not find a shorter periodicity at the W31230 data. If W31230 is indeed similar to SS 433, 
the lack of the short periodicity could be caused by a smaller inclination angle in our case. 
It would be therefore of special interest to obtain a spectrum of this interesting sta r, to see 
if there is any resemblance to the special spectral features of SS 433 (|Margonlll984^ . 

FoUowup spectroscopy turned out another interesting candidate, W21067, to be a fore- 
ground Calactic chromospheric star. We therefore suggest that the interesting cases found in 
M33 by photometry should all be followed by spectroscopy. With multi-object spectrographs 
available today on a few telescopes, such followup observations would require a relatively 
small amount of telescope time. 
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